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Abstract 



Suppose qi(x), i = 1,2 are smooth functions on R 3 and Ui(x,t) the solutions of the initial 
value problem 

d 2 U l - AUt - q l (x)U l = S(x, t), (x, t) eK 3 xl 
Ui(x,t)=0, fort<0. 

Pick R, T so that < R < T and let C be the vertical cylinder {(a;, t) : \x\ = R, R < t < T}. We 
show that if (U\, U\ r ) — (U2, U2r) on C then q\ = q 2 on the annular region R < \x\ < (R + T)/2 
provided there is a 7 > 0, independent of r, so that 

f \A s (q 1 -q 2 )\ 2 dS x < 1 I \ qi -q 2 \ 2 dS x , Vr e [R, (R + T)/2\. 

J\x\=r J|a:|— r 

Here As is the spherical Laplacian on \x\ = r. 
Partially supported by an NSF Grant DMS-0907909 
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1 Introduction 



Our goal is the study of a formally determined inverse problem for a hyperbolic PDE. Consider an 
acoustic medium, occupying the region IR 3 , excited by an impulsive point source and the response 
of the medium is measured for a certain time period at receivers placed on a sphere surrounding 
the source. We study the question of recovering the acoustic property of the medium from this 
measurement. 

Let q{x) be a smooth function on M 3 and U (x, t) the solution of the initial value problem 

U tt - AU-q(x)U = »irS(x,t), (x,t)GR 3 xM, (1.1) 
17 = 0, t<0. (1.2) 

Using the progressing wave expansion one may show that 

U(x,t) = 2 6 ^~ + u{x,t)H(t-\x\), (1.3) 



where u(x,t) is the solution of the Goursat problem 



utt ~ A« - q{x)u = 0, (x, t) G M , t > \x\, (1.4) 

u(x,\x\) = / q(ax) da. (1-5) 
JO 

The well posedness of the above Goursat problem is proved in |9| and improved in [11], though 
the result is not optimal; [9] has suggestions for obtaining better results and we will address them 
elsewhere. For completeness we restate the well posedness result. 

Theorem 1.1 (See {9] and 11 ). Suppose p > 0, and q is a C 8 function on the ball \x\ < p; then 



1.4), (1-5) has a unique C solution on the double conical region {(x,t) 6R 3 xR : \x\ < p, \x\ < 
t < 2p — \x\}. Further, the C 2 norm of u, on this double conical region, approaches zero if the C s 
norm of q, on \x\ < p, approaches zero. Also, if q is smooth then so is u. 



Below P =4 Q wm mean that P < CQ for some constant C. Let S denote the unit sphere 
centered at the origin. For any < R < T, we define (see Figure [TJ the annular region 

A := {x € M 3 : R < \x\ < (R + T)/2}, 

the space-time cylinder 

C = {(x,t) £l 3 xR : \x\ = R, R < t < T}, 

and 

K := {(x,t) G R 3 x M : R < \x\ < (R + T)/2, \x\ < t < R + T - \x\}, 
a region bounded by C and two light cones. 
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Figure 1: Geometry of the problem 



In our model the source is at the origin, the receivers are on the sphere \x\ = R and the signals 
are measured up to time T. Hence we define the forward map 



and our goal is to study the injectivity and the inversion of F. From the domain of dependence 
property of solutions of hyperbolic PDEs, it is clear that F(q) is unaffected by changes in q in the 
region \x\ > (R + T)/2. Hence the best we can hope to do is recover q on the ball \x\ < (R + T)/2. 

If q is spherically symmetric then the problem reduces to an inverse problem for the one di- 
mensional wave equation. In this case, recovering q on the region R < \x\ < (R + T)/2, from 
F(q), is done by the downward continuation method or the layer stripping method - see [16] and 
other references there. However, even in the spherically symmetric case (i.e. the one dimensional 
case), recovering gon |x| < R, from F(q) is more difficult since the downward continuation scheme 
is not directly applicable. It is believed that uniqueness does not hold for this inverse problem if 
T < 3R though explicit examples have not been constructed. If T > 3R, the question of recovering 
q on | a; | < R from (u,u r )\c was resolved by connecting this problem to one where the downward 
continuation method is applicable - see [8] and the references there. So it seems that in the gen- 
eral q case, recovering q over the region |x| < R will be harder than recovering q over the region 



Our main result concerns the problem of recovering q on R < \x\ < (R + T)/2 from (u,u r )\c- 
The downward continuation method does not apply directly in higher space dimensions since the 
time-like Cauchy problem for hyperbolic PDEs is ill-posed in higher space dimensions. Further, an 
analysis of the linearized problem shows that there could be singularities in q in certain directions, 
that is points in the wave front set of q, so that a signal emanating from the origin is reflected by 
this singularity in q, and the reflected signal never reaches the sphere |x| = R where the receivers 
are located - see Figure [2} Hence there should not be any stability for this inverse problem, unless 



we restrict q to a class of functions where singularities in q of the above type are controlled. In [13] , 

Sacks and Symes adapted the downward continuation method to apply to a slightly different inverse 



F : g H> (u\c,u r \c) 



R<\x\< (R + T)/2. 
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problem, with an impulsive planar source 5(z — t), with data measured on the hypersurface z = 0, 
where x = (y, z) with y G M 2 and z G R. They proved uniqueness for the linearized inverse 
problem when the unknown coefficient was restricted to the class of functions whose derivatives 
in the y direction were controlled by derivatives in the z direction. Later Romanov showed the 
inversion methods for one dimensional problems could be used for the existence and reconstruction 
for the nonlinear version of the Sacks and Symes inverse problem provided q(y, z) lies in the class 
of functions which are analytic in y in a certain sense, that is strong restrictions are placed on the 
changes in q in the y direction - see [lO] for details. We apply the technique in [13] to the uniqueness 
question for the problem of recovering q on on R < \x\ < (R + T)/2 from (u, u r )\c', we will have to 
impose restrictions on the angular derivatives of q. 

reflected ray 
4 




Figure 2: Reflection by a singularity in q 

For any i g I 3 we define r = \x\ and for i^Owe define 9 = x/r G S; hence x = r9. Define the 
radial vector field d r = r~ l x • V and, for 1 < i < j < 3, the angular vector fields Ojj = Xidj — Xjdi. 

Definition 1.2. Given 7 > 0, we define Q^(R, T) to be the set of all C 2 functions q(x) on the ball 
\x\ < (R + T)/2 with 

\\p\\HHs r ) + \\d r p\\m(s r ) < 7 {\\p\\m(s r ) + R?>lb(s r )) Vr e [R, (R + T)/2) 

where p(x) = f\ x ^ q(ax/\x\) da and S r is the sphere \x\ = r. 



So if q is a smooth function on |x| < (R + T)/2 with HpUff^Sr) + l|f?rP||z, 2 (S' T -) nonzero for every 
r G [R, (R + T)/2] then q G Q 7 where 

l|p||// 2 (S r ) + \\drP\\Hi(S r ) 
7 = max -TT—r — — n . 

re[R,(R+T)/2] \\p\\m(Sr) + \\°rP\\L2(S r ) 

Noting that d r p = q, using Garding's inequality on a sphere^], one may show that q G Q 7 * for some 
7* > if there is a 7 > so that 

\\A S q\\v(Sr) < l\\4hz(S r ), Vre[R,(R + T)/2] 



1 The Euclidean version is (6.8) on page 66 of [4]. Using a partition of unity argument and the Euclidean version, 
one may show that ||?||h 2 (s,.) — C r \\ Asg|| L 2( ST ,) with CV bounded if r is in a closed interval not containing 0. 
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where A5 is the Laplacian on S r . In particular, if q is a finite linear combination of the spherical 
harmonics with coefficients dependent on r then q £ Q^(R,T) for some 7 > 0. 

In section [2] we prove the following injectivity result using the ideas in |13| . 

Theorem 1.3. Suppose < R < T and qi,q2 are C 8 functions on R 3 . If F(q\) = F(q2) and 
Qi — Q2 £ Qy(R,T) for some 7 > then q\ = q2 on R < \x\ < (R + T)/2. 



One may tackle the problem dealt with in Theorem 1.3 using Carleman estimates also and one 
obtains a result which is stronger in some aspects and weaker in others. Using Carleman estimates 
one can prove uniqueness under slightly less stringent conditions on g - one needs controls on the 
L 2 norms of only the first order angular derivatives of p in terms of the I? norm of p, instead of 



on the second order angular derivatives required in Theorem 1.3 However, the price one pays is 
that the 7 cannot be arbitrary but is determined by R, T; further R cannot be arbitrary, but must 
satisfy R > T/2 and uniqueness is proved only for the values of q in an annular region R < \x\ < R* 
for some R* < (R + T)/2. This work will appear elsewhere. 



From Theorem 1.3 we can easily derive the following interesting corollary. 



Corollary 1.4. Suppose < T and (ft, #2 are smooth functions on R 3 which vanish in a neighbor- 
hood of the origin. If u\ and ui agree to infinite order on the line {(x = 0,t) : < t < T} and 
Qi — Q2 £ ^7(0, T) for some 7 > 0, then q\ = 52 on \x\ < T. 



We give a short proof of the corollary. If q± = q2 = in some small neighborhood of the 
origin then the difference u = u\ — u<i satisfies the standard homogeneous wave equation in a 
semi-cylindrical region 



{(x, t) £ R 3 x R : \x\ <S, \x\ < t < T - \x\}, 



(1.6) 



for some 5 > 0. Now, from the hypothesis, we have u is zero to infinite order on the segment of the 
t axis consisting of < t < T. Then by Lebeau's unique continuation result in [5] we have u = 
in the semi- cylindrical region given in (1.6). Hence u and u r are zero on the cylinder 

{(x, t) G R 3 x R : \x\ = 5, 5<t<T-5}. 



The corollary follows from Theorem (1.3) if the R and T in Theorem 1.3 are taken to be 5 and 
T — 5 respectively. 



We also have a uniqueness result for the linearized version of the inverse problem considered in 



Theorem 1.3; the result is for a linearization about a radial background. 



Theorem 1.5. Suppose qb(r) is a function on [0,oo) so that (ft(|x|) is a smooth function on R 3 ; 
further suppose Ub(r,t) is the solution of (1-4); (1-5) when q(x) is replaced by qb(\x\). Let q{x) be 
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a smooth function on M 3 and u(x, t) the solution of the Goursat problem 



u tt - Ait - q b u 
u(x, \x\) = 



gut, 
l 

q(ax) da. 





t > \x\, 



(1.7) 
(1.8) 



If (u, u r )\c = then q = on the region R < \x\ < (R + T) /2. 



This theorem holds with less regular q b and q; what is needed is enough regularity so that the 
spherical harmonic expansions of q, q b and u b converge in the C 2 norm. 



We next focus on the problem of recovering q on the region \x\ < R from (u, u r )\c when T > 3R. 
The linearized problem about the q = background, consisting of recovering q from (u, u r )\c, where 
u(x,t) is the solution of the Goursat problem 

uu — = 0, t > \x\, 
u(x,\x\) = / q(ax) da. 



As observed by Romanov, since T > 3R, we may recover q from (u,u r )\c fairly quickly. In fact, 
from Kirchhoff's formula (see (2l) expressing the solution of the wave equation in terms of the 
Cauchy data on C, we have 

, ,x f u r (y,t+\x-y\) ( 'u(y,t + \x - y\) u t (y,t + \x - y\)\ (y - x) ■ y 
u(x,t) = I — — h ( 1- — 1 — — -j ) — — — -^—dby. 



\x - y\ 



\x - y\ 



\x - y\ 



\x - y\ 



J\y-x\=R 

for all (x,t) with \x\ < t < R - see Figure [3j In particular we can express u(x, \x\) in terms of 
(u, ut, u r )\c and hence we can recover q. 



|y|=R 




s=t-|x-y| 



Figure 3: Kirchhoff's Formula 



For the original nonlinear inverse problem we show a partial uniqueness and stability result 
when one of the q is small. 
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Theorem 1.6. Suppose < 3R < T , M > and qi, i = 1, 2 are C 8 functions on \x\ < (R + T)/2 
wif/i Halloo < Af- Let ui be the unique solution of (1-4), (1-5) with q replaced by qi; then there is a 

(1.9) 



constant 5 > depending only on R,T and M so that if \\q 2 \\oo < 5 then 



\x\<R 



q 2 \ dx 4 I \u\ 
'c 



u 2 \ +\V(ui-u 2 )\ +\{u\-u 2 ) t \ dS Xjt ; 



the constant in (1.9) depending only on R,T,M. 



A weaker form of this result, requiring that \\qi\\ < 5 also, was given in [7j; a result similar 
to this weaker version was also derived in [12]. Later it was observed in [6j, for a similar type 
of problem, that the above proofs go through without the extra assumption that \\qi\\ < 6. We 
give this short proof of Theorem 1.6, in section |4j However, the original nonlinear inverse problem 
remains unsolved. 



2 Proof of Theorem 1.3 



2.1 Preliminary observations 



We need the following observations in the proof. For the angular vector fields we have [fiy, d T ] = 0, 
and [n -,n fc ,] = if = {k,l} but [n^fiifc] = fijy. Also |V/| 2 = / r 2 + r- 2 Z i<j (n ij f) 2 and if 
we define O = Yli<j then A = d 2 + 2r~ 1 <9 r + r~ 2 17 and [fi^, A] = 0. Also, for any i ^ j, since 
&ijf = Xidjf — Xjdif = dj(xif) — di(xjf) and XjXi — X{Xj = 0, by the divergence theorem, for any 
< i?i < R 2 we have 

n ij fdx = 0. (2.1) 



Ri<\x\<R 2 



Applying (2.1) to the zeroth order homogeneous extension of a function / on S, we conclude that 



for C 1 functions /, g on S 



n ij fdS = o, 



L 



fttijgdS 



L 



gQijfdS. 



(2.2) 



For i = 1,2 let u\ be the solution of (1.4), (1.5) when q = qi. Define Vi(x,t) = rui(x,t), 
Pi(x) = r Jq 1 qi(ax) da = qt{a6) da. Define v = v\ — v 2 , q = q± — q 2 and p = p\ — p 2 . Then we 
have 



Vtt - V rr - -p2^ v ~ Ql v = Q v 2, t > \x\ 

v(x, \x\) = p(x). 



(2.3) 
(2.4) 
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We are given that (v,v r ) are zero on C and we have to show that q = on R < \x\ < (R + T)/2. 
Note that since v = on C, we have p(x) = v(x, \x\) = on \x\ = R and hence for |x| > R we 
have p(x) = j r R q(cr8) da and hence to prove the theorem it will be enough to show that p(x) = 
on R < \x\ < {R + T)/2. 



We will attempt to carry out the proof which works in the one dimensional case. The limitations 
of this method when applied to the three dimensional case force the restrictions on q in the statement 
of Theorem 1.3 In the one dimensional case the angular terms are missing from (2.3) so the roles 
of r, t are reversible and one has sideways energy estimates which allow us to estimate the H 1 norm 
of v on t = \x\ in terms of the norm of v, v r on r = R and the I? norm of the RHS of (2.3). The 
H 1 norm of v on t = \x\ dominates the 1? norm of q on A and the L 2 norm of the RHS of the ( 2.3 ) 
is dominated by T — R times the 1? norm of q on A. So if T — R is small enough we obtain q = 
on A; then one combines a unique continuation argument with a repeated application of the above 
to prove that q = on A no matter what the T. 



In the multidimensional case the above argument breaks down because of the angular Laplacian 
in (2.3); all other parts of the argument work as in the one dimensional case. To carry out the 
above procedure we will need two estimates. The first is a standard energy estimate for the wave 
equation and the second is an imitation of a sideways energy estimate for a one dimensional wave 
equation in r, t where the roles of r and t are reversed. 



2.2 Energy identities 

For each p S [R, (R + T)/2], define (see Figure [4| the sub-region 

K p := {(x, i) £ I 3 x i : R < \x\ < p, \x\ < t < R + T - \x\}, 

the annular region 

A p := {x G M 3 : R < \x\ < p}, 

the vertical cylinder 

C p ■= {(x,t) £R 3 xl: \x\= p, p<t<R + T-r}, 

and for any function w(x,t) let w and w be the the restrictions of w to the lower and upper 
characteristic cones, that is 

w(x) = w(x,\x\), w(x) = w(x, R + T — \x\). 



We derive some relations which lead to the estimates we need. These relations are either the 
standard energy identity or a sideways version of it. Suppose w(x,t) satisfies 

Wit — w rr = F{x-, t), (x, t) G K. 



S 



T 



r+t=R+T 




O R A p (R+T)/2 



Figure 4: Sideways energy estimates 



Define the "sideways" energy (we will assume a sum over 1 < i < j < 3) 

J{p) ■= f r~ 2 (w 2 + w 2 + |V^| 2 ) dS X)t = f r- 2 (w 2 + w 2 + w 2 + r - 2 (n ijW ) 2 ) dS x>t 



R+T-p 
i p JS 

Multiplying the identity 



(w'f + w z r + w z + r~ 2 (VLijwf){pe, t) dOdt. 



2w r (wu — w rr —r 2 £lw — w) — 4r 2 £lijW r QijW + 2r 3 (flijw) 2 



l j vj r ij L l j UJ T ^1 \z L>iy I 

2 i „,.2 i m — 2/o ~..\2 i „..2\ i o/„,. „„ \ oo t 2„ 



-(^+w r +r (O^k;) + w ) r + 2{w r wt)t — 20ij(r w r flijiv) (2.5) 



by r , integrating over the region K p , using (2.2) and Stokes's theorem on a region in the r, t 
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plane, we obtain 

r~ 2 (2w r (F -w) — Ar~ 2 VLijW r VLijW + 2r~ 3 (Qijw) 2 ) dxdt 



K 



p r R+T-r 



-{wj + w 2 + r 2 (Qijw) 2 + w 2 ) r + 2(w r w t )t dt dr d6 



ISJR Jr 
(-T 



S JR 



(wf + w; + r- 2 (n ijW ) 2 + w 2 )(R6, t) dt d6 



R+T-p 



SJp 



S JR 



S JR 



(w 2 + w 2 . + r~ 2 (nijw) 2 + w 2 )(p6,t) dtdO 



(w 2 + w 2 r + r- 2 (n ijW ) 2 + w 2 - 2w r w t )(r6, R + T — r)drd6 
(w 2 + w 2 r + r~ 2 (Q, i: jw) 2 + w 2 + 2w r w t )(r9, r) dr dO 



J(R)-J(p)- [ r- 2 (w 2 + r- 2 (n ij w) 2 + id 2 )(x)dx 
JA P 

r~ 2 {w 2 r + r~ 2 {Slijw) 2 + w 2 )(x) dx. 



Hence 



J(p)+ f r~ 2 (\Vw\ 2 + w 2 )(x)dx + ! r- 2 (|Vu)| 2 + ?D 2 )(x)^+ t 2r~ 5 (n 
= J(R) + / r~ 2 [2ww r + 4r~ 2 QijwMijW - 2Fw r ) dxdt, R< p < 

JK 



w) dx dt 



R + T 



This is the sideways energy identity we need, 
t 



r+t=R+T 




r+t=R+T 




(R+T)/2 



(2.6) 



Figure 5: Standard energy estimate 

Next we derive the standard energy identity for the wave equation. For any s £ [R, T], define 
(see Figure [5]) the domain 

K s = K(~){(x,t) £ M 3 x M : R < t < s}, 
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H s the horizontal disk obtained by intersecting K with the plane t = s, that is 

H s = KC\{t = s}, 

whose projection onto the plane t = is the annular region 

A s := {x e M 3 : R < \x\ < min(s, R + T - s)}. 

Next, we define the "energy at time s" for every s E [R, T] - the definition depends on s < (R+T)/2 
or not because the geometry changes - see Figure [5] For s G [R, (R + T)/2], we define (summation 
over 1 < k < I < 3) 

E(s) := / r~ 2 (w 2 + w 2 + w 2 + r~ 2 (0 M w) 2 )(x, s) dx 
and for s€[{R + T) /2, T] we define 

E(s):= r~ 2 (w 2 + w 2 + \Vw\ 2 )(x, s) dx + r- 2 (w(x) 2 + \(Vw)(x)\ 2 ) dx. 

JA s J R+T-s<\x\<(B,+T)/2 

First take s < (R + T)/2; multiplying the identity 
2u; t (u;i t -u; rr -r~ 2 ilu; + t(;) = (w 2 + w 2 + w 2 + r~ 2 (f} feZ u;) 2 ) t - 2(w t w r ) r - 2r~ 2 9 Jk i(w t Slkiw) (2.7) 
by r~ 2 , integrating over the region K s , and using (2.1), we obtain 

/ 2r~ 2 w t {F + w)dxdt 

= E(s) + 2 / r~ 2 (w t w r )(x,s)dS x dt 

= E(s) + 2 f [ r- 2 (w t w r )(x,s)dS x dt- [ r - 2 (w(x) 2 + \(Vw)(x)\ 2 ) dx. 

JR J\x\=R JA s 



Next take s € + T)/2,T]; multiplying (2.7) by r , integrating over the region K s , using (2.1) 



11 



we obtain 



/ 2r 2 w t {F + w)dxdt 

= / r~ 2 (w 2 + wj + \Vw\ 2 ) dx 

+ / r~ 2 (w 2 + w 2 + w 2 - 2w t w r + r~ 2 (n k iw) 2 )(x,R + T - \x\) dx 

J R+T-s<\x\<{R+T)/2 

- [ r~ 2 (w 2 + w 2 + 2w t w r + r- 2 (Q H w) 2 ){x, \x\)dx + 2 [ [ 

J A JR J\a 

= / r~ 2 (w 2 + w 2 + \Vw\ 2 )dx+ / r~ 2 (ii){x) 2 + \(V'5))(x)\ 2 )dx 

JH° JR+T-s<\x\<(R+T)/2 



r' 2 {w{x) 2 + \{Vw){x)\ 2 )dx + 2 I I r- 2 (w t w r )(x,t) dS x dt 

A JR J\x\=R 

= E{s)- [ r- 2 {w{x) 2 + \(Vw)(x)\ 2 )dx + 2 [ [ r~ 2 (w t w r ){x,t) dS x dt. 

J A JR J\x\=R 

Hence, in either case, that is for any s £ [R, T], we have 

E(s) < / r~ 2 (w{x) 2 + \{Vw){x)\ 2 )dx + 2 / r~ 2 w t (F + w) dx dt + / r~ 2 {w 2 t + w 2 ) dS x dt. 
J A JK" Jc 

(2.8) 

2.3 Uniqueness 



r (wtw r )(x, s) dS x dt 

\x\=R 



We now show that if v and v r are zero on C then q = on A. We apply (2.6) to v = v\ — v% note 
that F = vtt — v rr — r~ = qiv + qv2 and J(R) = because the Cauchy data of v is zero on C. 
Hence 

J(p)+ / r~ 2 (v 2 + |Vu| 2 ) < / r~ 2 (vv r +4r~ 2 QijV r QijV -2v r (q!V + qv 2 ) 
JA P JK P 

4 [ r- 2 (v 2 +v 2 .+r- 2 (n ij v) 2 + q 2 + r- 2 (n ij v r ) 2 ) 
JK P 

P p / rR+T—r \ /■ 

J(r)dr+ / r- 2 q 2 (x)( / dt) dx + r~ 4 (%u r ) 2 



< f P J{r)dr + (T -R) [ r~ 2 q 2 (x)dx+ [ r~ 4 (^ r 
JR J A JK 



i2 



with the constant associated to ^ being ci = 4max(l, ||<Zi||.l°°(a)) H u 2||i,«>(in)> Hence, by Gronwall's 
inequality 

J(p)+ f r - 2 (\Vp\ 2 +p 2 )^(T-R) [ r- 2 q 2 (x)dx+ [ r-\n ijVr ) 2 , R < p < 

JA P J A JK 1 

(2.9) 



12 



with the constant being C2 = c\e Cl( - T R \ In particular 

J(p)4(T-R) [ r~ 2 q 2 (x)dx + [ r~ 4 (QyW r ) 2 , R 
J a Jk 



<P< 



R + T 



and taking p = (R + T)/ 2 in (2.9) we have 



r~*(\Vpr +P ) =4 (T-R) / r~ l q 2 (x)dx + / r~ 4 (^-u r )' 
a J a Jk 



(2.10) 



(2.11) 



with the constant c-i- Integrating (2.10) w.r.t p over [R, (R + T)/2] we obtain 



/ r- 2 (v 2 + v 2 + \Vv\ 2 ) 4 (T-R) 2 [ r~ 2 q 2 (x) dx + (T — R) [ r- 4 (fiy« r ) a . (2.12) 
Jk Ja Jk 



So we can combine (2.11), (2.12) into 



f r- 2 (v 2 + v 2 + \Vv\ 2 ) + [ r- 2 (p 2 + \Vp\ 2 ) 4 (T - R) [ r' 2 q 2 + [ r^^Vrf 
Jk Ja Ja Jk 

with the constant being C3 = (1 + T — R)c2- 



(2.13) 



The equation (2.13) would have been enough to prove Theorem 1 in the one dimensional case, 
because |Vp| 2 > p 2 = q 2 and the last term in (2.13) would not be there. Then by taking T — R 



small enough we could have absorbed the second term on the RHS of (2.13) into the LHS and we 



would have proved the theorem for T close to R. Then a unique continuation argument would 
prove the theorem for all T > R. However, in the three dimensional case we do have the last term 



in (2.13) which cannot be absorbed in the LHS because it involves second order derivatives of v 



we will estimate it in terms of p using the standard energy estimate for the wave operator. 



with Q. Note that from (2.3) and (2.4) we have 



Fix an i,j pair with i < j. We apply (2.8) to the function w = noting that fij,- commutes 



with 
and 



Wtt - Wrr - -^2^ w = F 



F(x,t) := qiw + (Qijqi)v + (flyg)«2 + Q^ijV2- 



(2.14) 



w(x,\x\) = (n ijP )(x). (2.15) 
Further, since the Cauchy data of v is zero on C, so the Cauchy data of w is zero on C. Hence 



from (2.8) we have 

E(s) < [ r- 2 ((n ijP ) 2 + \Vn ijP \ 2 ) + I r~ 2 (w 2 + w 2 + F 2 ) 
J a Jk* 

4 [ r- 2 ((n ijP ) 2 + \Vn ijP \ 2 )+ [ r- 2 (w 2 + w 2 + v 2 + q 2 + (n ljq f 

J A JK" 

4 [ E(t)dt+ [ r- 2 (p 2 + \Vp\ 2 + \Vn ijP \ 2 )+ [ r~ 2 v 2 
Jr J a Jk 
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with the constant being C4 = 2max(l, (R + T) 2 , \\qx\\oo, W^ijllWoo, Halloo)- So from Gronwall's 
inequality we have 



r~V, R< s<T 



E(s) 4 [ r- 2 (p 2 + |Vp| 2 + \Vn ijP \ 2 ) + / 
J A JK 

with the constant being C5 = c^e C4 ^ T ^ R \ Integrating this w.r.t s over the interval [R, T] we obtain 
f r- 2 (w 2 + w 2 t + \Vw\ 2 ) < c 5 (T - R) ( [ r - 2 (p 2 + |Vp| 2 + |V%p| 2 ) + f 

JK \J A JK 



-2 2 

r v 



hence, since w = QijV, 



f^r-\n ijVr ) 2 <c 5 R- 2 {T-R) ^ r - 2 (p 2 + \Vp\ 2 + \Vn ijP \ 2 ) + J^r~ 2 v 2 ^j . (2.16) 



Using this in (2.13), we have 



r- 2 v 2 + / r- 2 (p 2 + \Vp\ 2 ) 
K J A 



4(T-R) [ r- 2 (p 2 + \Vp\ 2 + \Vn ijP \ 2 ) + (T-R) I r~V (2.17) 

J A JK 



with the constant cq = max(c3, C3C5R 2 ). However, q is in Q 7 so 



f (R+T)/2 r 

r- 2 \V(n ijP )(x)\ 2 dx= / / (Vtt ijP )(r6) 2 d6 dr 

a jr J\e\=i 

f(R+T)/2 f - 

<7 / r 2 (p 2 + \Vp\ 2 )(re) dO dr 

jr J\e\=i 



<^(R + T) 2 [ r- 2 {p 2 + \Vp\ 2 ). 
J A 



Using this in ( |2.17[ ), we see that p = on A if T — R is small enough - depending on 7, cq and 

f R 

Jo 



R + T. Now v(x, \x\) = p(x) and v = on \x\ = R so p = on \x\ = R, that is q(a9) da = for 



all unit vectors 6. Hence 

j q(a8) da = 0, R < r <T 
Jr 

which implies q{x) = when R < \x\ < T, provided T — R is small enough. 

Actually, adjusting the height of the downward pointing cone, what we have shown is the 
following: there is a 5 > dependent only on j,R,T, || || c" 1 (A) ? 1 1 ^2 He 1 CRT)? so that ^> f° r some 
R* £ [R, {R + T)/ 2], v and v r are zero on the cylinder 

{(x,t) : \x\=R*, R* <t<R* + 26}, 

then q = on R* < \x\ < R*+5, with the obvious modification in the assertion if R*+5 > (R+T) /2. 
We use this observation to prove that q = for any R, T. 
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Since v and v r are zero on C, then from the above claim, we have q = on R < \x\ < R + 5. 
Let u = u\ — U2 where ui,uz are solutions to (1.4), (1.5) for q = q±,q2- Then, u satisfies the 
homogeneous equation 

utt — Au — q\u = 

over the region K p where p = R + 5. Now u and u r are zero on C, and q\ is independent of t, so by 
the Robbiano-Tataru unique continuation theorem (see Theorem 3.16 in [KKL01]) we have u = 
in the region K p ; in particular u and u r are zero on C p and hence v, v r are zero on C p . Now repeat 
the above argument, except R is replaced by R + 5; this argument repeated will complete the proof 
of Theorem 11.31 



3 Proof of Theorem 1.5 



Let {(/> n (x)}^ =1 be a sequence of homogeneous harmonic polynomials on M 3 so that their restrictions 
to the unit sphere S form an orthonormal basis on L 2 {S) - see Chapter 4 of 15 . Let k(n) be the 
degree of homogeneity of (fi n . Then q(x) and u(x,t) have spherical harmonic decompositions in 
L 2 (S) given by 



r k{n) q n {r 



q(r9) = qn(r)r k ^M9), u(r8, t) = J] u n (r, 1)^^(9) 
q(rO) MO) dO, 



n=l n=l 

where 

r k ^u n (r,t)= [ u(re,t)<j> n (6)de. 
J\9\=x 

Since u and q are smooth, we may s1iok|^] that q n (r) and u n (r, t) decay as n~ p for large n for any 
positive integer p, uniformly in r, t. Hence the series also converge in the C 2 norm. 



To prove the theorem, it will be enough to prove that q n (r) = on R < r < (R + T)/2 for all 
n > 1. One may show that for sufficiently regular / (see page 1235 of \\\) 



A [f(r,t)r 



fc(n), 



.(*) 



Mn) Me)(ftt-frr 



2k(n) 



'fr 



hence, using (1.7), (1.8), the u n (r,t) are solutions of the one dimensional Goursat problems 
,2„. «2.. 2fc(n)-2. 



d t u n 



9~ u n ■ 



u n (r, \r\) 



drUn 
1 



qbu n = q n Ub, 



t > r 



a k ^qn(ar)da. 



o 



The hypothesis of the theorem implies that u n (R,t) and (d r u n )(R,t) are zero for R < t < T. So 
repeating the standard argument for one dimensional hyperbolic inverse problems with reflection 



data, as in 17 , or repeating just the sideways energy argument in the proof of Theorem 1.3 without 



the complication of the angular terms, one may show that q n (r) = for R < r < (R + T)/2. 



2 Use the definition of q n and u n , observe that the <f) n {9) are eigenvalues of the spherical Laplacian, and use the 
Divergence Theorem on 5* to transfer the Laplacian from the cj> n to q or u - see Theorems 2 and 4 in [14| . 
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4 Proof of Theorem 



1.6 



t=T 




Figure 6: Transmission data problem 

Let (see Figure |6| B denote the origin centered ball of radius R in M 3 , D the region 
D := {(x,t) e K 3 x R : |x| < R, \x\ < t < T}, 
and as before C the cylinder 

C := {(x,t) : \x\ = R, R<t< T}. 



Let Ui, i = 1, 2 be the solutions of ( |1.4[ ) , (|L5j) when q = q%, define q = q\ — q% and u = u\ — 112- 
Then u satisfies 



u a - Au- qiu = qu 2 , (x, t) G D 

u(x,\x\) = / q(ax) da. 
Jo 



(4.1) 
(4.2) 



Then, restricting attention to the cylindrical region B x [R, T] , from [3] we have the following 
stability estimate for the time- like Cauchy problem (note T > 3R): there is a constant C\ dependent 
only on M, R, T so that 



K-'OIIhMB) + K(->*)lli2(B) ^ °l {\\9U2\\h(Bx[R,T]) + IMI/P(C) 



+ Wl3(oJ. r < t <t. 

(4.3) 

Next, if we multiply (4.1 ) by ut and use the techniques for standard energy estimates (backward in 



time) on the region |x| < t < R, we obtain 

\u{x)\ 2 + |Vu(x)| 2 dx < C 2 I JJ^ ^ ^ \qu 2 \ 2 dxdt + |K, R)f H i {B) + KO, R)\\h {B) j (4.4) 



where n(x) = u(x, \x\) and C2 depends only on M,R. Hence, combining (4.3), (4.4) we obtain 



J \ u ( x )\ 2 + |Vu(x)| 2 dx < C 3 (||<?u 2 || 2 2(D) + \\u\\ 2 Hl(c) + |Kllia(0)) ( 45 ) 



16 



where C3 depends only on R, T, M. Now ru(x) = J Q r q(s8) ds, hence q(x) = (ru) r = u + ru r . So 
q 2 <2(u 2 + r 2 u 2 ) < 2 max(l, R 2 ){u 2 + u 2 ) < 2 max(l, R 2 )(u 2 + | Vu| 2 ), 

and 

Ikllia(B) ^ C 4 (||g«2|lia(fl) + \\u\\hi( C ) + [l"rllia(C7)) ( 4j 
with C4 dependent only on R,T,M. Finally, using Theorem we have 

II<^2||l 2 (D) < ||«2|U<»(D) lkllL 2 (D) < A/"(T, 1 1 Q-2 1 1 00 ) lkl|L 2 (D) 

where the Halloo norm is over the region \x\ < (R + T)/2. Since M(T, Halloo) goes to zero as H^H 
approaches 0, we can choose a 5 > so that 

cwV(t,|Moo) < J 



if 1 1 <Z2 1 1 00 < 5; note that this <5 will depend only on R, T, M. Using this in (4.6), we conclude that if 
1 1 <?2 1 1 00 < 5 then 



(4.7) 

with C5 dependent only on R, T, M. 
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